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CT -1 Abstract 
Ph 

<T^ I derive sharp semiclassical asymptotics of J |e/,(x, y, 0)| 2 oj(x, y) c/x dy where e/,(x, y, r) 

ch is the Schwartz kernel of the spectral projector of Magnetic Schrodinger operator and 

uj(x, y) is singular as x = y. I also consider asymptotics of more general expressions. 

S 

Introduction 

> 

This paper is a continuation of [Ivrl2] and I consider 

o 



(0.1) / = J J u(x, y)e(x, y, r)^ 2 (x)e(y, x, r)^i(y) dx dy 



where e(x,y,r) is the Schwartz kernel of the spectral projector E(r) of the (magnetic) 



Schrodinger operator 



S3 (0.2) A= Uj2Pjg jk (x)P k -v), Pj = hD j - f ,Vj 

in IR d and r = 0, h — > +0 while /i — > +oo. 
Further, 

(0.3) w(x, y) = Q(x, y; x - y) where function ft is smooth in 6(0, 1) x 6(0, 1) x B(R d \ 0) 
and homogeneous of degree — k (0 < k < d) with respect to its third argument ^ 



1 ' In other words it is Michlin-Calderon-Zygmund kernel. 
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and smooth cut-off functions ipi, fa- 
It follows from [Ivrl2] that 1=1+ 0(h 1 ~ d ~ K ) as fj, = 0(1) with X defined by the same 
formula but with e(x, y, r) replaced by 

(0.4) e y w (x, y , r) = f (2nh)~ d [ e ^H-y.O 

J g{y^)<v{y)+2r 

Then the standard rescaling technique implies the same asymptotics but with the remainder 
estimate 0(fih 1 ~ d ~ K ) provided 1 < /i = o(/? _1 ). 

Let d = 2. Then in the general case it is the best remainder estimate possible while 
0( / u/?~ 1 ) is the best possible remainder estimate for 

(0.5) J= J e(x, x, r)^(x) dx; 

one needs to consider constant magnetic field 

(0.6) F = g-Hd Xl V 2 - d X2 V{), g = detig^y 1 

and g Jk = const, V = const. In this case spectrum consists of Landau levels 

(0.7) l -{{2n + l)^hf - V), neZ + , 

of infinite multiplicity and e(x,x, r) = e MW (x, r), 

(0.8) e MW (x, r) = i- ^ 9(2r + V - {2n + l^h^^fy/g 

n>0 

is a Magnetic Weyl expression; see [Ivrl]. 

On the other hand, in the generic case the remainder estimate for (0.5) is o(fi~ 1 h~ 1 ) 
and the principal part is J e MW (x, r)^(x) dx if F does not vanish and /i < h -1 ; otherwise 
the remainder estimate for 0.5 is o(fi~ 1 ^ 2 h~ 1 ); moreover in the latter one can consider 
h^ 1 < fi < h~ 2 and the remainder estimate would be the same but the principal part would 
be 0(/i -1 /? -3 ) and the formula could be more complicated ([Ivr7, Ivr8]). 

Now my purpose is to get the sharper remainder estimate for (0.1) under the same 
conditions. This is a very daunting task since for (0.5) periodic trajectories were the main 
source of trouble and they were broken in the generic case; for (0.1) loops are also the source 
of trouble, and in the generic case former periodic trajectories generate a lot of loops (see 
figures 1 and 2). 
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As the result, while an asymptotics with the sharp remainder estimate but with the 
principal part given by very implicit Tauberian formula (1.4) is a rather easy corollary 
of [Ivrl], Chapter 6 (done in Section 1), the deriving of an asymptotics with the sharp 
remainder estimate and rather explicit principal part is much more difficult. 

In Section 2 I consider the weak magnetic field case when /i is not very large and 
therefore magnetic drift is relatively fast and the next winging is distinguishable in the 
quantum sense from the previous one and replace T in the Tauberian expression by e/i -1 
with certain error estimate; combining with [Ivrl2] I get a remainder estimate (sharp as 
/i < (h\ log/?!) -1 / 6 but not very shabby for /x < (h\ log /?|) -1 / 3 ) with non-magnetic principal 
part. 

Section 3 is devoted to the strong magnetic field case /i > h^ 1 / 2 ^ 5 when the reasonable 
remainder estimate could be derived by the method successive approximation and the 
unperturbed operator is a model operator, admitting explicit calculations. 

Finally, in Section 4 I consider the strong magnetic field case when both approaches are 
combined. 

The results of sections 2-4 are not always sharp or very very explicit, but could be made 
either sharp or completely explicit. But some calculations are left to the readers. 

1 Estimates 

1.1 Tauberian Formula 

So I am considering operator (0.2) where g-* k , Vj, V are smooth real-valued functions of 
x G M. d and (g^) is positive-definite matrix, < h <C 1 is a Planck parameter and fih < 1 
is a coupling parameter. I assume that A is self-adjoint operator. Then simple rescaling 
x i— > fix, h i— > fih, /x i — > 1 leads us to the remainder estimate 0(fih 1 ~ d ~ K ) which is completely 
sufficient for applications to Multiparticle Quantum Theory. 

However I want to improve this remainder estimate under generic non- degeneracy con- 
dition. I consider only the most sensitive case d = 2 rather than d = 3. 

1.1.1 In this paper I assume that 



(1.1) 
(1.2) 
(1.3) 



V>e, 
|F| > e, 
\VV/F\ > e 



where F is an intensity of magnetic field. Then 



1 ESTIMATES 



4 



Proposition 1.1. Under conditions (1.1) — (1-3) the contribution of zone {\x — y\ > C^}, 
to the remainder is 0(/j,~ 1 h~ lr y~ K ) while the main part is given by the same expression (0.1) 
with e(x,y, 0) replaced by its standard implicit Tauberian approximation with Ixe/i 

(1.4) e T (x,y,0) d = f h' 1 f F t ^ T {x T {t)u{x,y,t)) dr. 

J — oo 

Proof. Here and below U(t) = e' h 1{A is the propagator of A and u(x,y, t) is its Schwartz' 
kernel. 

Consider expression (0.1) with uj(x,y) replaced by c^ 7 (x,y) which is a cut-off of c<;(x,y) 
in the zone {|x — y| x 7} and with the original ipi-ipi replaced by 1. Let us replace one 
copy of e(x, y, r) by e(x, y, r, r') = (e(x, y, r) — e(x, y, r')) with r' < r and the second copy 
by e(x, y, t") and denote the resulting expression by / 7 (r, r', r"). 

Now let us use decomposition 

(1.5) w 7 (x, y) = 7 ~ d ~ K y ^i, 7 (x, z)^ 2 , 7 (y, z) oz 

where d = 2 now. 

Then / 7 (t, t', t") does not exceed 

(1-6) ^C 7 -1|^(r,r> y || 1 

y 

where E(t,t') = E(t) — E(r'), tpj are real-valued 7-admissible functions supported in C07- 
vicinities of z,- and B(zj, 2Q7) are covering of our domain of multiplicity not exceeding Co. 
Here I used that ||E(r")|| = 1. Since E(r, r') is a positive operator and ipj = ipj, one can 
replace trace norm by the trace itself and get 

(1.7) ^C 7 - K Tr^E(r,r> i = C^Tr E(r, r> 

j 

with $ = Ej^ 2 - 

Further, I know from the standard theory [Ivrl] that under conditions (1.1)-(1.3) 

(1.8) \\E{t,t')^\\i < Ch- 2 (\T-T'\ + CT^h) Vr.r'e [-e,e], T = efi. 
and therefore 

(1.9) mr,r',r")\ < C^h' 2 ^ - r'\ + CT^h) 

in the same frames and therefore due to the standard Tauberian arguments I conclude that 
the contribution of zone {|x — y| x 7} to the Tauberian remainder estimate does not exceed 
C/x _1 /?~ 1 7~ K which implies the statement immediately. □ 
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However I need to consider zone {|x — y| < C7}, complementary to one above. Assume 
that 

(1.10) ft K (z) = D zP K -u + fi«-i,o 

j 

with index at ft showing the degree of the singularity. Then 

(1.11) lv k (x,x- y)^ 7 (x -y) = Y^ Px,(w«-lj^ 7 ) + W «^7 + 

j 

where ?/> 7 = ip((x — y)7 _1 ) with ip supported in 6(0, 1) and equal 1 in 6(0, \) while 
is defined similarly with ip' supported in 6(0, 1) \ 6(0, |) and the last term gains 1 in the 
regularity. 

After integration by parts expression / K?7 , defined by (0.1) with ft replaced by ft^, 
becomes 

(1.12) -^Ysjj ^-u(^y)(hD XJ )(e(x,y,r) ■ e(y,x,r)) dxdy 

j 

plus two other terms: the term defined by (0.1) with kernel Q' K j of the same singularity 
k, without factor h' 1 and supported in the zone {|x — y| > 7/2} and the term defined by 
(0.1) with kernel ft^_i 7 - without factor h^ 1 and of singularity k — 1. 

The former could be considered as before yielding to the same remainder estimate 
0(fi~ 1 h 1 ~ d 'y~ K ). To the latter I can apply the same trick again and again raising power 
(and these terms are treated in the same manner (but simpler) as I deal below with (1.12). 

So, one needs to consider (1.12) and thus, denoting the second copy of e(y,x, r) by 
f(y,x, t) and without using that they are equal 

(1.13) (hD Xj )(e(x,y,r)-f(y,x,r)) = 
(hD Xj e(x,y,T))f(y,x,r) - e(x, y, r)(f(y, x, r)(hD Xj Y) = 
(P lx e(x,y,r))f(y,x,r) - e(x, y, r)(f(y, x, r)P/ x ). 

I remind Pj = hDj — fiVj(x) and Pj = —hDj — fiVj(x) is the dual operator. I also remind 
that if e(x,y,r) and f(y,x, r) are Schwartz kernels of E(r) and F(r), then P,- X e(x, y, r) 
and f(y,x,T)Pj x are those of PjE(r) and F(t)Pj. 

Therefore I am interested in the expressions of the type 



1.14) h 1 JJ u K ^ l {x,y)e{x,y,T)f{x,y,T)ip 1 dxdy. 
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If k < 1 then replacing e(x, y, r) (and f(y, x, r)) by its standard Tauberian expressions one 
gets an error not exceeding Ch^ 1 x fi~ 1 h~ 2 'y 1 ~ K because ||P,-E(r)|| < Co, ||P//~(t)|| < Co 
where F(r) is an operator with the Schwartz kernel f(x,y, r) and also because 

(1.15) Y^\\<PjPjE{ry)<pA\i <J2yjPMr,T')\\ 2 -\\E(r,T') n \\ 2 < 

j j 

j j 
J2^VjPjE(r,r')P^j + J2^VjE(r,T') V j 

j j 
< Ch- d (\T-r'\ + Cn~ x h) Vr,r'e[-e,e] 

which also easily follows from [Ivrl]. 

So, in this case one gets remainder estimate 0(fi~ 1 h 1 ~ dr y~ fi + fi~ 1 h~ dr y 1 ~ K '} which is 
optimized to 0(/i -1 /? 1 ~ d ~ K ) as r x h. 

On the other hand, as 1 < k < 2 one can apply the same trick again since I did not 
use the fact that e(., ., .) and f(., ., .) coincide; then I arrive to the same estimates with Pj 

replaced by PjPk or even by P J d = Pj\P } \ ■ ■ ■ Pj,: 

(1.16) Tr P J E(t,t')(P j )* < Ch- 2 (\r-T'\+ fi-H) Vr, t' G [— e, e]. 

Finally, note that as k ^ 1, (1.10) is always possible. Further, as k = 1 decomposition 
(1.10) is possible as well provided one adds term x(x)|x — y| _1 with an appropriate coeffi- 
cient. On the other hand if k = 1 and u(x, y) = x(x)|x — y| _1 then (1.10) is also possible 
but with u j(x,y) = x(x)(x / - — y/)|x — y| _1 log |x — y|. 

So I arrive to 

Proposition 1.2. Let conditions (1.1) — (1-3) be fulfilled. Then 

(i) As < k < 2 and either k ^ 1 or k = 1 and u is replaced by uo — x(x)|x — y| _1 with an 
appropriate coefficient x(x), with the error 0(^ 1 h 1 ~ d ~ K ) one can replace e(x,y, r) by its 
standard Tauberian expression (1.4) in the formula (0.1) for I . 

(ii) As k = 1 and to = x(x)|x — y| _1 , with the error 0(fi~ 1 h 1 ~ d ~ K \ \ogh\) one can replace 
e(x,y,r) by its standard Tauberian expression (1.4) in the formula (0.1) for I. 

Remark 1.3. The arguments above show that in an appropriate sense one can consider 
arbitrary K6l and even in C 

(ii) Can one prove the similar result for l m defined by (0.6), [Ivrl2] with m > 3? 
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1.1.2 Since one needs only (1.16) rather than (1.1)-(1.3) and (1.16) holds as (1.3) is 
replaced by a weaker non-degeneracy condition (see [Ivr9]) 

(1.17) V/F has only non- degenerate critical points 

1 arrive to 

Proposition 1.4. Let conditions (1.1), (1.2) and (1.17) be fulfilled. Then (1.4) and state- 
ments (i),(H) of proposition 1.2 hold. 

Remark 1.5. Under certain assumptions (see [Ivr6]) this result could be generalized for 
d > 4. However in calculations I will need conditions (1.1)-(1.3). 

2 Calculations: Weak Magnetic Field 

I am going to assume from now on that d = 2, // < h^ 1 and conditions (1.1)- (1.3) are 
fulfilled. In this subsection I assume that the magnetic field is weak enough. I remind 
that classical particles move along cyclotrons which are circles of the radius x /i -1 (and 
respectively ellipses if g jk = const) as g-* k = d~jk, V = const, F = const but which drift in 
more general assumption with the velocity /i~ 1 V( V/F)-*-. This illustrates the difficulty I 
am facing: trajectories 2 -* are coming back. In this subsection I consider weak magnetic field 
approach when one gets sharp remainder estimate even if one ignores returning trajectories. 

2.1 Isotropic Approach 

2.1.1 I want to replace Xr(t) by Xr'(t) with T' = efi -1 in the Tauberian formula and 
to estimate the corresponding error; in the correct framework this would be equivalent to 
using non-magnetic Weyl approximation for e(x, y, 0) (see (0.3) [Ivrl2]). 

So let us consider classical Hamiltonian dynamics at Figure 1. Let us notice that 
while the trajectory is not periodic due to non-degeneracy condition, it is self-intersecting. 
However, generic point on the trajectory is not a point of the self-intersection, and on each 
trajectory winging the number of self-intersection points is x /j and the length of the loop 
is at least x ^r 1 (actually typical loop is of the length x 1 but I will make a more precise 
statement and use it later). 

The trajectory as F = const and V is linear is hypercycloid as on figure la; in more 
general case the trajectory drifts along some curve and the cyclotron radius changes as on 



2 ' Even if one calls them trajectories they are projections of actual trajectories. 



2 CALCULATIONS: WEAK MAGNETIC FIELD 



8 




(a) hypercycloid (b) perturbed hypercycloid 



Figure 1: Classical trajectories 

figure lb but still one can consider x on the first winging and direct drift at this moment 
along Xi. 

Let us parametrize x on the first winging by € [0 , 2tt] ; without any loss of the generality 
one can assume that <fi £ [0,7r/2]; one can reach it by some combination of reflections, 
reverting time direction and permutation of x and y. So, x belongs to the bold arc on 
figure 2: 




(a) NP is the North Pole, = (b) Loop analysis: right side 



Figure 2: Classical trajectories: NP is on the top of the first winging 

So, let us consider some fixed trajectory and point x(0) on it with < <fi < 7r/2; this 
fixes also £(0) 3 ). 

3 ' In conformal coordinates, where (g^) is proportional to Euclidean metrics, cf> is an angle between 
V(V/F) and ( Pl ,p 2 ). 
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Note first that in virtue of [fvrl] 

Proposition 2.1. Let <t» t be drift flow on the energy level (see [Ivrl]). Let ipi be supported 
in Cqh~ 1 -vicinity ofx. Then 

(2.1) \F t ^h-wX,4t-*)( 1 -Mx))<x,y,t)il>i{y)\ < Ch 5 Vr : |r| < e^ 1 
as ipi = 1 in Co/i™ 1 -vicinity of<t> to (x). 

I claim that 

Proposition 2.2. Let Q = Q(x, hD) be operator with (e, fie)- admissible symbol where 

(2.2) Cifi'Hl log h\f' 2 < e < e 

Then for \t\ < Co, Q t = U(—t)QU(t) is also (e, fie) -admissible operator with supp Q t = 
U' t (supp Q) with the corresponding Hamiltonian flow ty t . 

Proof, (i) Let us rescale x i— > x/i, h i— > h~ = fih, fi \— > 1, e \— > log h\yi 2 . 

As F = 1 the proof is really easy since then e 2nlh toA is a standard h-FlO corresponding 
to symplectomorphism W 27r ^-i different by 0(/i _1 ) from identical and thus 

(2.3) e 2*;ft-M = g/M- 1 ^"^ 

where /. is /i-PDO commuting with /4; then 

(2.4) e''^ lt4 = e '-ft- 1 fi e ft"»-M 

with t' = hn, n = [fit /{In) + 1/2J (and then |t'| < CyU 2 /? < 1 and t" = fit - 2irn (and then 
< c) and both of these operators are standard ft-FIOs. 

It is a bit more complicated in the general case but one can always assume that F(x) = 1 
and then (2.3) still holds and then the same arguments hold as well. □ 



2.1.2 Let us calculate the distance from x on the upper-right quarter of 0-th winging to 
the nearest point y on the right part of n-th winging (n e Z \ 0). One can see easily that 

(2.5) The distance from x to the nearest point y on the right part of n-th winging is 
r(0, n) x fi~ 2 \n\ sin provided | sin <p\ > C yU _1 |n|. Otherwise £((f), n) x fi~ 3 n 3 . 



2 CALCULATIONS: WEAK MAGNETIC FIELD 



10 



This condition means that the distance x pr 1 sin 2 to the horizontal line passing through 
the North Pole (0 = 0) is larger than the distance one measures; meanwhile the deviation 
of the enveloping curve from this horizontal line is 0(p~ 4 n 2 ) which under this condition is 
much smaller. As n = 1 such distance is shown on figure 2b above by blue arrow. 

Therefore I conclude that point (x, £) is distinguishable from each point W t (x, £) with 
e/i" 1 < \ t\ < C , residing on the right-halves of the windings provided 

(2.6) e/i~ 2 max(0, p" 1 ) > e = C(p~ l h\ log h\) 1/2 

where I took the smallest possible e. 

Note that one can satisfy (2.6) for any provided 

(2.7) p<e{h\\ogh\Y^; 

let us assume temporarily that this is the case. 
Let us introduce (e, /xe)-admissible partition 

(2.8) J2 Q J = I 
and define 

(2-9) q ± =J2 q j< Q± (p)= E q j 

jeJ± jeJ±{ P ) 

where refers to elements residing on the right-halves of the trajectories and -/ ± (p) refers 
to those elements of J^ 1 which are in p-vicinity (with respect to 0) of the poles; finally 
J(p) with p = eo/i -1 refers to "polar caps" (p-vicinities of the poles) and I put them to 
arbitrarily. 

So, Tauberian formula (1.4) for ej becomes 

(2.10) e T (x,y,0) = h- 1 ^ / F t ^ T (x T {t)Q ]x u(x, y, t)Ql y ) dr 

U,k)eJxj J -°° 

and therefore 

(2.11) f T {x, y) = e T {x, y, 0) - e T/2 (x, y, 0) = 




T- 1 Ft^h-w(xT{t)Q Jx u{x,y, t)C?y | r=0 

(j.k)eJxJ 
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with x{t) = — x(2t), x(t) = — /'t _1 x(t); here J = J x J contains all the pairings but 
later I also consider J = J ± x J^ 1 and J = ^ x J T . 
Let us consider 

(2.12) /J- = y w(x, y)/ f T (x, y)e(y, x, 0) c/xt/y. 

Due to the same arguments as before 

(2.13) If one replaces in (2.12) uo by its cut-off in the zone {|x — y\ > 7} then the result 
would not exceed CT~ 1 h~ 1/ ~f~ K . 

We apply this estimate with 7 x p~ l T. Note that as T e [Q),e/i] zone {|x — y| < 7} 
with 7 = eoP~ 1 T provides a negligible contribution into (2.12) due to proposition 2.1 and 
therefore 

(2.14) |/^| < CJ-^^h- 1 as T G [C ,e/i]. 

On the other hand, if I replace in (2.12) uo by its cut-off in the zone {|x — y\ < e 1 p~ 1 T 2 } 
and pick J = x I get a negligible result as well due to proposition 2.2 and (2.5). 

So, let us consider 7 e [ei/i -1 7~ 2 , eyU -1 7~], 7" e [e//~\ Co] and let us replace in (2.12) 
uo by its cut-off in the zone {|x — y\ x 7}. Then one can also replace by ^(p) with 
p = cpT~ lr y; the error will be negligible again due to proposition 2.2 and (2.5). 

To estimate the resulting expression I need an inequality 

(2.15) \\Q{p)E{r,r')Q{p)\\ 1 < Cph-\\r - r'\+ p-'h) 

which I will prove a bit later. Due to (2.15) and our standard analysis such modified 
expression (2.12) with J = J* 1 x does not exceed CT~ 1 p^~ K h~ 1 (with 7 = p~ l Tp) i.e. 

(2.16) CT- 1 - K p 1 - K p K h- 1 . 
Summation with respect to p from eo T to 1 results in 

(2.17) C7- 1 - K (7" 1 - K + 1 + 6«i| log T\)p K h- 1 x C{T~ 2k + T' 1 -^! + b Kl \ log T^p^" 1 

and summation with respect to T from 7" to 7~i results in the same expression (2.17) with 
T replaced by 7" . In particular, for 7~o = e^" 1 I get 

(2.18) C(p 3K + p 2K+1 + p 3 b Kl \ \0gpDh- 1 

So far I replaced in just one copy of e-r(x,y, 0) in (0.1) T = ep by the smaller value (and 
only in J ± x ^ pairs. However exactly the same arguments work for the second copy as 
well. 

So, I arrive to 
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Proposition 2.3. Let conditions (1.1) — (1-3) and (2.7) be fulfilled. Let us replace in the 
Tauberian formula for e-r(x,y, 0) (plugged into {0.1)) T = e /i by T 6 [e^ -1 , e//]. JTien 

(i) As 7~o > Co i/je error does not exceed the right-hand expression of (2.14); 

(ii) As T < Cq the contribution to the error of all pairs (J, k) e J + x J + U J~ x J~ does 
not exceed (2.17). 

(iii) In particular, as T = e^T 1 this contribution does not exceed (2.18). 
Proof of (2.15). Proof is standard based on the standard calculation of 
(2.19) F^ h - lT (x T (t)Tr{Q(p)U(t))); 

I leave details to the reader. □ 

Actually one can draw conclusions as (2.7) is violated but our present arguments are 
too crude anyway. 

2.1.3 Now let us consider J = J* x J T ; due to (2.13) one needs to consider T < Q only; 
so the error in pairs would not exceed (2.17) with T = 1 and one should not bother 

to get better estimate anyway. 

Let analyze the left halves of trajectories (Figure 3). 




(a) Right part (b) Left part 



Figure 3: Classical trajectories: loop analysis: left side 

In the arguments of the previous subsubsection condition (J, k) 6 J*- x J*- was used 
only to dismiss certain distances between x and y as impossible. In particular the same 
estimates imply that 

(2.20) If one replaces in (2.12) u by its cut-off in the zone {|x — y| > ep~ 2 } the result would 
not exceed CT~ l n 2K h . 
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2.1.4 In this subsubsection I consider the case eo < (fi n < n/2 — eo, assuming that e < ji 2 
i.e. 

(2.21) /i<e(/?|log/?|)- 1/3 . 

Let <j) n be the location of n-th intersection (or intersection of original quarter with n-th 
winging). Consider point x which is closer to n-th intersection than to any other intersection 
point. Consider y on m-th winging. 

Let us analyze case n ^ m first. Then unless \<f>„ — 4> m \ < ei, the distance between x 
and y is x ^r 1 and one should not be concerned since such points are covered by (2.20). 
On the other hand, if \(f> n — m | < e 1 , the distance between x and y (as n ^ m) is at least 
fi~ 2 and such pairs are covered by (2.20) again. 

So, one needs to consider only m = n; now one needs to consider contributions of the 
pairs (J, k) connected by a trajectory, and "gravitating" to the same intersection point. 
Let us joint all the elements "gravitating" to n-th intersection point and residing on the 
distance not exceeding 7 from it with 7 e [Ce, c/i -2 ]; more precisely let us define Q^ and 
Q nr/ as the corresponding sums of Qj. Then by the standard methods of [Ivrl] one can 
prove easily that 

(2.22) ||(? n , 7 f(r,r')<?„, 7 ||i< C^h' 2 (|r - r'\ + fi^h) . 

Furthermore, if either supp Qj or supp Q k was on the distance x 7 from the intersection 
point, then the distance between supp Qj and supp Qk is also x 7 and therefore the con- 
tribution of all such pairs to the error does not exceed C/xn~ 1 7 1 " K . Then summation with 
respect to 7 from e to Cjj,~ 2 results in 

(2.23) C^ 2K - 2 + e x ~ K + 5 Kl | log^el)^" 1 
as T x 1. Meanwhile 

(2.24) The contribution to the error of pairs when both elements Qj and Qk are supported 
in Ce-vicinity of the intersection point does not exceed Cfieh^ 1 ^ K . 

The proof of (2.24) repeats arguments leading to proposition 1.2 with (fi framing Q^E(r, t')Q^ 
from both sides. Clearly 0/ and Q n do not commute well (scales are incompatible) but one 
does not need a commutation here. 

Together (2.23) and (2.24) imply that the contribution of one "tick" (number n) to the 
error does not exceed 



(2.25) 
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However there are x p ticks and one needs to multiply by p resulting in expression 

(2.26) Cfi^h- 1 + C^h- 1 ' 2 -^ log/)] 1 / 2 . 

The first term here is exactly as in (2.20) with T x 1. Therefore I arrive 

(2.27) If one replaces in the Tauberian expression in pairs (J, k) £ x J T with at least 
one element, residing in zone {min /gZ \<p — nl/2\ > e }, T = C by T = e with arbitrarily 
small constant e, it would cause the error in (0.1) not exceeding (2.26). 

On the other hand, since in the zone in question the distance between connected elements 
is x as \t\ < e we conclude that replacing T = e by some smaller value T would cause 
the error in (0.1) not exceeding CT' 1 p K h~ 1 . Combining this with (2.27) we arrive to 

(2.28) If one replaces in the Tauberian expression in pairs (J, k) £ J ± x J T with at least 
one element, residing in zone {m\n i e i\(f) — n I /2\ > e }, T = Co by some smaller value T, it 
would cause an error in (0.1) not exceeding 

(2.29) CT^^h- 1 + Cp^h' 1 + Cfi 3/2 h' 1/2 ~ K \ log h^ 2 ; 
in particular, as T x /i -1 one gets 

(2.30) Cfi^h- 1 + C^ 2 h-V 2 - K \ log h^ 2 . 

2.1.5 Let us consider now <p n x p < e (figure 3a). More precisely, let us consider 
contribution of pairs (J , k) e x J T such that both elements Qj and Qk are in p- vicinity 
of the pole and at least one of them is on the distance p from the pole. As I want to use 
Tauberian formula with T < e p I need three jumps: from T = Co to T = Cop, from 
T = Cop to T = top and from T = ep to desired T. As T x p there would be only two 
jumps and as p > CqP there would be only one jump. 

One can see easily that in the first jump the distance is at least e Q p^T 2 and applying 
our standard arguments I conclude that the error does not exceed 

(2.31) C7- 1 p(A*" 1 T 2 )- K /7- 1 x CT-^p^h- 1 ; 

similarly in the third jump the distance is at least e p~ 1 p 2 and the error does not exceed 

(2.32) CT- 1 p(p- l p 2 )- K h- 1 x CJ- 1 p 1 - 2K p K h- 1 . 



The second jump as one could see is more tricky. 
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Note that in zone in question <p n x p~ l n and \<p n — m | x p~ l \n — m\ as m ^ n. Since 
branches intersect under angle 2<p n I conclude that as n ^ m the distance is measured 
approximately along the vertical and it is x p~ 2 \m — n\p i.e. at least p~ 2 p. 

Repeating the same arguments as before I conclude that the contribution to the error 
does not exceed 

(2.33) CT- 1 p{p' 2 p)- K h- 1 x CT- K p 2K h-\ 

As m = n the distance again is measured approximately along the vertical and it is 
x 7p. Therefore x and y are distinguishable as p7 > Ce i.e. 7 > 7' with 

(2.34) i = Cpep-\ 

Again I need to assume that 7 < ep~ 2 , i.e. that p~ 2 > 7' = Cpep -1 and it is the case even 
for p x /i -1 under assumption (2.7). 

So, the contribution of such pairs is estimated by CT -1 p r ){p r y)' K h~ l and summation 
with respect to 7 from ep^ 1 to p~ 2 results in (2.23)-like expression 

(2.35) CT' 1 p- K [p lK ' 2 + e 1 - R p K - 1 + b Kl \ log p^ep-^ph' 1 . 

Meanwhile the contribution of pairs with both elements in 7'-vicinity of the intersection 
point does not exceed 

(2.36) C{pep- 1 )T- l h- l - K x Cp~ 2 p 1 / 2 ^ 1 ' 2 ^ . 

So contribution of one "tick" to the error does not exceed (2.35) + (2.36) but there are x pp 
ticks and taking in account that p x T we arrive to (2.26)-like expression 

(2.37) Cp-^p^h- 1 + Cp-y /2 /T 1/2 "1 log h] 1 / 2 . 

Finally, summation with respect to p > T results in (2.37) with p replaced by T 

(2.38) CT-^p^-H- 1 + CT-V 3/2 /7" 1/2 "1 log /)| 1/2 
while summation of (2.31) with p > 7" results in 

(2.39) CT~ 2K p K h^ + CT^p^h- 1 + Cr-V^T 1 ! log r|6 Kl/2 . 

Note that (2.31) with p = T is just a first term in (2.39). 

In particular, plugging into (2.39) T x p' 1 and adding to (2.38) I get 

(2.40) Cp^h" 1 + Cp 5/2 h~ 1/2 ~ K \ log h\ 1/2 + Cp^h' 1 + C/i 2 ^ 1 ] log T|6 Kl/2 . 
So, I conclude that 
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(2.41) If one replaces in the Tauberian expression in pairs (J, k) E J ± x J T with at least 
one element residing in zone {min/ e ^ \<f> — 7r/| < eo} T = Co by some smaller value T, it 
would cause an error in (0.1) not exceeding (2.38) + (2.39); in particular, as T x pr y one 
get (2.40). 

2.1.6 Let us consider now zone n/2 — e < 0„ < n/2 (figure 3b). Then there exists n (the 
last intersection) such that 

(2.42) (p- n < C fi~ 1/2 , ip„ x fi~ 1/2 {n - n) 1/2 as n ^ n, £ n = \ip„ - ip n+x \ x 1 (p~ \ 

with ip n d = (tt/2 - <f) n ); 

so ticks are actually longer than before. 

Since windings intersect under angle 2<p n I conclude that as n ^ m the distance is 
measured approximately along the horizontal and it is x pr 2 \m — n\ sin((/9„ + y? m ). Then x 
and y are distinguishable as long as it is larger than e. This is always the case as /i~ 5 / 2 > e 
i.e. 

(2.43) /i< epilog/)!)- 1 / 4 . 

Note that one needs to consider only n, m with ip n , <p m e [0, e ] (since the other pairs are 
already covered). Then the contribution to an error does not exceed 

C f jL- 1 T- 1 h- 1 (jji- 2 ( p m )- K x Cfjt^h-^n-ml-*' 2 x C^-H' 1 -^ 2 x C^h' 1 

m<n m<n 

since k < 2. This is exactly the first term in (2.26). 

Consider now m = n < n. Then the distance again is measured along the horizontal and 
it is x /i -1 |y? — (p n \<£n (<f = ( 7r /2 — 0)- Then x and y are distinguishable as /i _1 |y? — <f n \<fn > 
Ce. 

So, the contribution of such points is estimated by another (2.23)-like expression 

(2.44) C^T-H- 1 J 7- K c/ 7 x C^h-'^-Hr' + ie^Y^ + d^oghl) 

where integral is taken from 7 = eip~ 1 to 7 = fi~ 1 £ n . 

Meanwhile the contribution to the error of pairs with both x and y residing in zone 
{\<f — (p n \ < tp^e} does not exceed CfiT~ 1 h~ 1 ~ K e(f~ 1 ; one can prove it easily by the same 
methods as before. Adding to (2.44) I get the contribution of pairs residing near n-th tick 

(2.45) Cfi 2 ^- 1 ^ 1 + C^eh- 1 -^- 1 x CjV^/r 1 + /le/T 1 -^ x ^ 1/2 (n - n)- 1 ' 2 
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where I rewrote the first terms in (2.44) due to (2.42); other terms of (2.44) are dominated 
by the last term in the left side of (2.44). Note that the first factor on the right is exactly 
(2.26) while the second factor sums to C with respect to n. 
Finally, let us consider m = n = n or, more precisely, zone 

{\<p\ < e fi~ 1/2 }. Then the 

above arguments still work as long as <p > C(fie) 1 ^ 2 = C(/ih\ log/?!) 1 / 4 and I arrive to the 
term not exceeding C / u 2k_1 / 2 /?~ 1 plus contribution of zone {\<p\ < C(fih\ log/?]) 1 / 4 }, which 

does not exceed 

(2.46) epilog/?!) 1 / 4 /?- 1 -" 

which in turn does not exceed Cfi^ 1 h^ 1 ^ K under condition (2.7) and the second term in 
(2.26) otherwise. 

Therefore I arrive to 

(2.47) Claims (2.27), (2.28) (with an extra term Cfi^ 1 h^ 1 ^ K in estimates) are valid in the 
"near equator" zone as well. 

2.1.7 So, when one replaces in (1.4) T = Cq by T = e// -1 , the error in u would not 
not exceed (2.18) + (2.40) + C/j l ~ 1 h~ 1 ~ K (since (2.26) + (2.30) is lesser); removing terms 
dominated by others I arrive to 

(2.48) C(/i 2k+1 + Qi 3 *)/! -1 + C/i 5/2 /T 1/2 - K + Cfi-'h- 1 -*. 
Thus I arrive to 

Proposition 2.4. Let conditions (1.1) — (1-3) be fulfilled. Then under assumption (2.7) 
/ is given by the standard Tauberian formula with T = e/i -1 with an error not exceeding 
(2.48). 

This result is very crude: condition (2.7) is too restrictive and estimate (2.48) is not 
sharp enough. In the next subsection improving this approach I get better results. 

2.2 Anisotropic Approach 

In this subsection I will improve results of subsubsections 2.1.2 and 2.1.5 to weaken condition 
(2.7) and to improve the error estimate; I will analyze even the case of condition (2.21) 
violated. 
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2.2.1 The improvement is based on anisotropic partition elements. More precisely I will 
take ^-scale along trajectories and e-scale in the perpendicular direction (with respect to x, 
and matching partition with respect to £) with 

(2.49) e<^<e fi-\ e>C ^ 2 , 

(2.50) eq > C^ l h\ \ogh\ 

where the second condition is intended to counter the curvature of trajectories and (2.48) 
will be logarithmic uncertainty principle. 

However, one needs to overcome certain obstacles: first, one needs to explain what it 
actually means (how to quantize such symbol). Second, scales of different boxes are not 
compatible and they are not also compatible with the original (2.8) partition. Finally, such 
operators do not work nicely with FIOs and one needs to overcome this obstacle. 

To quantize a symbol supported in a straight box is easy: it is just a standard Weyl 
quantization. In the general case one can transform any rotated box into straight one by 
rotation which is a metaplectic operator; then ne can nicely quantize it by Weyl. On the 
other hand, exactly the same result would be achieved if one quantized the original symbol 4 -* 

One can see easily that such operators nicely commute with operators corresponding 
to (q, /x^)-boxes; moreover two such operators corresponding to two rotated boxes with the 
angle 9 between corresponding axis also nicely commute if 

(2.51) e 2 > C\s\n9\ ■ \ogh\. 
In particular, 

(2.52) If Qi, Q 2 are operators with the symbol supported in a rotated box 81,82 and if 
doubled boxes do not intersect, then Q 1 Q 2 = 0. 

So, I need to study propagation of singularities: basically I need to prove that 

(2.53) Let \t\ < T < Co. Let (y,fj) = ^(x.^), be a quantization of the symbol 
supported in the rotated (e, q, jiq, fie)- vicinity of (3c, £) and Q 2 be quantization of the symbol 
equal to 1 in C (e, q, fiq, /ie)-vicinity of {y,rf). Then (1 — Q2)U(—t)Q 1 = 0. 

This is not a theorem (yet) because conditions to e and q except (2.49) and (2.50) are 
missing. On the figure 4a the corresponding vicinities are a filled and an empty rectangles 
and here may be much more than one winging. 

4 ) Weyl quantization of symbols obtained by linear symplectic transformations in the phase space co- 
incides with the metaplectic transformation (corresponding to this symplectomorphism) of the quantized 
symbol, in particular with (x,£) (qx, g^ 1 ^) which I apply here; means a transposed matrix. 
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(a) Anisotropic propagation (b) to 2.2.2 



Figure 4: Anisotropic boxes 

To prove (2.53) I need to elaborate the nature of operator L in (2.4), (2.5); I remind 
that I rescaled x i— > x/i, h \— > h = fih, /i \— > 1, but did not rescale respectively t — > t/x, 
however I used \t'\ < T. 

Assume for a moment that g-* k = const, and Vj, V are linear. In this case 7t x e' f Hi (x, £) = 
<t>t'(x) where <lv(x) = x+vt' and pT l v is the speed of the drift. Then e 27r ' M 4 = where 
72^1/ means a shift u(x) i— > u(x — 2nv) (modulo a scalar factor). 

Then obviously the general case e' h t>L = T t i v !F t i where T t > is ft-FIO with the Hamilto- 
nian map / + 0(/i~ 1 7" 2 ). So, evolution is described by a product of metaplectic operator 
corresponding to rotation and generated by a quadratic Hamiltonian (i.e. Hamiltonian of 
the form q(x, fihD) where q is quadratic form on R 4 ), operator corresponding to the linear 
phase shift and generated by a linear Hamiltonian (i.e. Hamiltonian of the form £(x, fihD) 
where £ is linear form on IR 4 ) and the operator corresponding to the symplectomorphism 
which is almost /. 

One needs to investigate how such operator acts on the operator which is fi-quantization 
o the symbol supported in (e' , , , e')-box; here I made a replacement q i— > q' = qp, 
e e' = efi; then (2. 49), (2. 50) become 

(2.54) e' < q' < e , e' > C q' 2 , e'q' > Ch\ \ogh\. 

Now rescaling again x 1 \— > xi/p, x 2 i— > x 2 , D 1 \— > Di, D 2 i— > x 2 /p, ft i— > ft/p with p = e/q I get 
the last symplectic map / + 0(q) and an isotropic settings. Then I arrive to 
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Proposition 2.5. Let e' > C{n~ x h\ log h\fl 2 i.e. 
(2.55) e > Cp^h 1 ' 2 

and let Q[ be a quantization of the symbol supported in the rotated (e',q',q',e') -vicinity of 
(x, £) and Q' 2 be a quantization of the symbol supported in the rotated Cq(e' , q', q f , e')-vicinity 
of the same point (x, £). 
Then (/ - Q^t'Q 1 ! = 0. 

I would not need more analysis of propagation for the sake of arguments of the next 
subsubsection (see figure 4b) since as |x — y| > q one can apply analysis of subsubsection 
2.1.2 without any modification. 

However, for analysis of x J T pairs in later I will need a bit more analysis since the 
boxes will be tilted rather than almost parallel. 



2.2.2 I replace (2.8), (2.9) by different partitions while considering the different partitions 
but Q + and Q still have sense. In this subsubsection we consider x pairs. Then I 
should take 

(2.56) e = eofi" 1 Tp, q = ^{Tp) 1 ' 2 

where p x T + | sin 0| as before, the first equality is needed since the distance is x p~ x T p 
and the second follows from the second inequality in (2.49) 5 - ) ; then due to (2.50) one needs 
to assume that 

(2.57) pT> C(ph\\ogh\) 2 '\ 

Then (2.54) is fulfilled as well. 

I remind that p > T and therefore (2.57) holds everywhere as 

(2.58) T > C{ph\\og h\f' 3 . 
On the other hand (2.57) holds for p x 1 as 

(2.59) T > C{ph\\og h\) 2 ' 3 . 
In particular, (2.58), (2.59) are fulfilled with T = ep' 1 as 

(2.60) M < log /?|) 1/4 



5) 



Obviously one should take the largest values possible. 
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and 



(2.61) 



V<e(h\\ogh\)- 2 / 5 



respectively. 

Since in the calculations of subsubsection 2.1.2 the actual size of of e did not matter, 
only the sheer fact of being disjoint did, I conclude that 

Proposition 2.6. Let conditions (1.1) — (1-3) be fulfilled. Let us replace in the Tauberian 
formula for e T (x,y, 0) (plugged into {0.1)) T = e p by T e [e/i~\ Co]. Assume that (2.56) 
holds. Then 

(i) The contribution to the error of all pairs (J, k) G J + x J + U J x J does not exceed 



(ii) In particular, as T = e/i 1 this contribution does not exceed (2.18). 

On the other hand, if only (2.59) holds but (2.58) fails this analysis remains valid (only) 
in zone 



and then I arrive to 

Proposition 2.7. Let conditions (1.1) — (1-3) be fulfilled. Let us replace in the Tauberian 
formula for ej{x,y, 0) (plugged into (0.1),) T = e p by T e [e/i" 1 , Co]. Assume that (2.57) 
holds but (2.56) fails. Then 

(i) The contribution to the error of all pairs (J, k) e J + x J + U J x J with p > p T 
does not exceed 



(2.63) CT- l -*(p\-* + 1 + 5 Kl | log TD^h- 1 x 

C(T- 2 {ph\ log /,|)2(i -)/3 + r -i-«(! + 5k1 | | 0g T\))p K h- 1 



(2.17). 



(2.62) 



p>p T ^CT-\ph\\ogh\yH 



(2.64) 




Finally, in this latter case I need to add contribution 



(2.65) 



Cpp T h 1 K 



CT-\p 5 / 2 h\\ogh\) 2 l z h 



, — 1 — K 



of the polar caps {p < pj to the error and I arrive to 
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Proposition 2.8. In frames of proposition 2.7 

(i) The contribution to the error of all pairs (J, k) e J + x J + U J x J - does not exceed 

(2.66) cr- 1 ^ 5 / 2 ^ log /)|) 2 / 3 /?- 1 - K + CT~ 1 - K p K h- 1 

(ii) In particular, as T = e/i" 1 this contribution does not exceed 

(2.67) C{p 4 h\ log hi) 2 / 3 /)- 1 -" + C/S^h- 1 

2.2.3 Consider now contribution of x J T pairs. In this case the size of q does matter 
from purely geometrical point of view. The other problem is that the boxes introduced in 
subsubsection 2.2.1 are mutually tilted now. 

To avoid this problem let us partition Q + and Q~ in different manner. I need to satisfy 
(2.49), (2.50). I also need to satisfy e > qp. Therefore I select 

(2.68) e = C{p~ l h\ log /)|) 1/2 p 1/2 , ? = Cifi-Hl log h^p' 1 ' 2 . 

Here I analyze only case p < eo since results of subsubsections 2.1.4 and 2.1.6 could not be 
improved by considering of anisotropic boxes. 
Now 

(2.69) x and y are distinguishable as 7 > q. 
Really, one can prove easily 

(2.70) Let conditions (2.49), (2.50) and e > qp be fulfilled. Then symbol bo9 with b(x, f ) = 
(3(e~ 1 x 1 , q~ 1 x 2 , s' 1 ^), regular symbol (3 and 9 = 9 p a rotation of both x and £ by angle 
p is also quantizable 6 ) and the corresponding metaplectic transformation of the resulting 
operator coincides with the quantization of b. 

Then encapsulating both rotated and mutually tilted non-intersecting boxes into two 
C (q,£,£,q) boxes (2.69). 

Note that q < p" 2 (which is the length of the tick) as 

(2.71) P>Pi = Cp 3 h\\ogh\. 

Note that (2.69) is satisfied for all p (I remind that p > T and thus in the final run p > /i -1 
under (2.60). On the other hand (2.69) holds for p x 1 under (2.21). 

6 ' I remind that I use [i~ x h quantization. 
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Then as in (2. 35), (2. 36) the contribution of one tick satisfying (2.71) to the error when 
one replaces T = Cop by T = eop does not exceed 

(2.72) CT^h' 1 [p~ > 2K + /i^/r K ) x Cp^/r 1 (p-y^ 1 + {ph\ log h\) 1/2 p- 1/2 h~ K ^ 

and contribution of x pp ticks does not exceed 

(2.73) Cp/T 1 (VV*- 1 + {»h\ log h\) 1/2 p' 1/2 h- K ^ 

and summation with respect to p > 7" returns the same expression with p replaced by 7" 

(2.74) Cph- l [r- K p 2K ~ l + (p/)| iog/71) 1 / 2 r^ 1 / 2 /?—). 

Noting that contribution of {t : \t\ > T, t ^ p} (as it was derived in subsubsection 2.1.5) 
does not exceed 

(2.75) Cp K T^h- 1 (1 + T x - 2k + 5 K i | log T\) 
I arrive to 

Proposition 2.9. Lei us replace in the Tauberian expression in pairs (J, k) £ x J T with 
at least one element residing in zone {min/ e z \4> ~ < e o} value T = Co by some smaller 
value T satisfying 

(2.76) T > max(ep-\ Cp 3 h\ log h\). 

Then this would cause an error in (0.1) not exceeding (2.74) + (2.75). 

In particular, as T x p~ y and (2.60) is fulfilled the error does not exceed 

(2.77) Cp^/r 1 + Cp 2 h^ 2 ^ K \ log h\^ 2 + Cp^/T 1 + Cp^h' 1 ^. 

Here (2.75) with T = p" 1 brought only one term which is not necessarily dominated 
by (2.74) with T = p" 1 or Cp^ 1 h~ 1 ~ K . Therefore in comparison with (2.40) I weakened 
condition (2.7) to (2.60) and gained factor p~ 1//2 in the second term. 

Assume now that 

(2.78) ep' 1 < T<p 1 = Cmin(Cp 3 /?|log/?|,l). 

Again let us replace first T = C$p by T = e p. Then the contribution of zone {| sin 0| x p} 
with e T < p < pi does not exceed Cp" 1 x ph~ x ~ K x h~ l ~ R ; I remind that the polar caps I 
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already covered in subsubsections 2.2.1, 2.2.2 and excepted them from analysis here. Then 
summation with respect to p results in 

(2.79) Or^llogT/pil. 

Meanwhile, as p < (h\ log Z?!)" 1 / 3 contribution of zone {| sin 0| > pi} does not exceed (2.74) 
with T x pi which one can see easily that it is less than (2.79). On the other hand, as 
P > {h\ log /?|)~ 1//3 one should take pi = 1 and this latter zone disappears. 

Further, as p~ x p\t\ < h\ \ogh\ one should reconsider contribution of zone {\t\ ^ p, p < 
C(ph\ log/?!) 1 / 2 }. This however falls inside of the "polar cap". So, I arrive immediately to 

Proposition 2.10. Let us replace in the Tauberian expression in pairs (J, k) e x J T 
with at least one element residing in zone {min /gZ |</> — irl\ < e } (with the exception of polar 
caps) value T = Co by some smaller value T satisfying (2.78). Then this would cause an 
error in (0.1) not exceeding (2.78) + (2.75. 

In particular, as T x p~ x and (2.60) fails this error does not exceed 

(2.80) Ch" 1 "*] \og{p 4 h\ log h\)\ + Cp K+1 h-\ 

2.2.4 Now I want to combine results of two previous subsubsections as T = ep~ l . I remind 
that each of them contains two statements, as T exceeds some critical value or below it; 
these critical values are Tj* = C(ph\ log/?!) 1 / 3 and 7" 2 * = Cp 3 h\ \ogh\. Let us compare them 
and also with ep^ 1 . 

Proposition 2.11. Let conditions (0.3) and (1.1) — (1-3) be fulfilled. Let us replace in the 
Tauberian formula T = Cq by T = ep^ 1 . 

(i) Assume that condition (2.60) holds. Then the error does not exceed (2.18) + (2.74) + 
(2.75); in particular, as T x p~ l an error does not exceed 

(2.81) Cp 2K+1 h- x + Cp^h' 1 + Cp 2 h' 1/2 ' K \ log h\ 1/2 + Cp^h' 1 ^. 

(ii) Assume that but (2.60) fails but (2.61) holds. Then 

(a) as T > 7" 2 * = Cp 3 h\ log ft] an error does not exceed (2.18) + (2.74) + (2.75), 

(b) as T*>T> 7"* = C(ph\ log/?!) 1 / 3 an error does not exceed (2.18) + (2.79) + (2.75), 

(c) as T{ > Tep^ 1 an error does not exceed (2.66) + (2.79) + (2.75). In particular, as 
T x p^ 1 an error does not exceed (2.67). 
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Proof, (i) Note that T x * > T| iff (2.60) holds. However in this case 7~| < e/i _1 which means 
that the error does not exceed (2.18) + (2.74) + (2.75). As T = efi' 1 reducing dominated 
terms I arrive to (2.81). 

(ii) Assume that (2.60) fails but (2.61) holds. Then T* < e^ 1 < 7? < 7" 2 * and 

(a) As (2.21) holds and T G [7" 2 *, C ] an error does not exceed (2.18) + (2.74) + (2.75); 

(b) As T e [T*, T*} an error does not exceed (2.18) + (2.79) + (2.75); 

(c) As T e [e/i-\ 77] (2.66) + (2.79) + (2.75). As T = e/i" 1 reducing dominated terms I 
arrive to (2.67). 

□ 



2.3 Main Theorem 

Now I can prove the main theorem of this section: 

Theorem 2.12. Let conditions (0.3) and (1.1) — (1-3) be fulfilled. Then I is given by the 
standard non-magnetic Weyl expression 7 ^ with the remainder estimate given by (2.81) if 
(2.60) holds and by (2.67) i/ (2.60) fails but (2.61) holds. 

Proof. In view of proposition 2.11 one should consider I defined by a standard Tauberian 
expression with T = e^T 1 . Scaling x 1— > fix, fi 1— ► 1 and h 1— > fih I arrive in view of the proof 
of Proposition 2.6 [Ivrl2] to expression 

(2.82) J ~ x mn h- 2 - K+m+n ii m \ 

m,n 

however since for real kernel u I must be real as well, and complex-conjugation is equivalent 
to fj> 1 — —fi, only even powers of /i are allowed; then modulo 0(fj, 2 h~ K ) I arrive to the same 
expression with \x = 0. Note that fi 2 h~ K is well below the announced remainder estimate. 

Due to [Ivrl2] again this expression equals to the standard Weyl expression modulo 
0(h~ K ). □ 



7 > See (2.27)-(2.29), (2.7) [Ivrl2]. 
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3 Calculations: Strong Magnetic Field 

3.1 Calculations for a Model Operator 
3.1.1 Consider first a model operator 

(3.1) A = ^([hD 2 ) 2 + {hD 1 - fix 2 ) 2 - W - 2wc 2 ), W = const > 0, v = const. 
Then 

(3.2) A = ^ ((/?D 2 ) 2 + {hD 1 - fix 2 + n^vfiih - (1/1/ + n~ 2 v 2 - fih) - Iv^hD^j 

and therefore with (2.3) holds with L = fiB, 

(3.3) B = n({W + n~ 2 v 2 - fih) - 2vn~ 1 hD l ) . 

Let us consider the corresponding spectral problem with 1/1/ = (so including a constant 
part 1 1/1/ of the potential into a spectral parameter): 

(3.4) l - [[hD 2 f + (/A - fix 2 ) 2 - 2vx 2 ) e = Ae 

with e = e(x, y, A). 

Making a (partial) unitary 8 ^ Fourier transform with respect to xi one gets 

(3.5) X - ([hD 2 ) 2 + (/»& - /ix 2 ) 2 - 2vx 2 ) e = Ae 
with e = e(fi, x 2 , 771, y 2 ) or equivalently 

(3.6) ^(/?D 2 ) 2 + (/if ! - ^x 2 + vfi- 1 ) 2 ^ e = (2A + 2i/ A T 1 /7f 1 + >/ 2 /i" 2 )e; 
therefore 

(3.7) e = n 1/2 h-V 2 T n (^ 2 h-y 2 (-^h^+x 2 -vn- 2 ))x 

fi:(2n+l)^<2A+2i/M~ 1 ^l + *'V~ 2 

r n (^h-V2 { _ fI -l h {: i+y2 _ vfT^e-^Sfa-Th) 

8 > So, with the factor {2k)- 1 ' 2 . 
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where T n are (real orthonormal) Hermite functions. Therefore excluding again |H/ from 
the spectral parameter one gets 



(3.8) e(x,y,0) = (2vrrV/ 2 /,- 1 /^ f 

„ J (In 



X 

(2n+l),uh< W+2v f j,- 1 h^ 1 + v 2 ^- 2 



T n (/X 1/2 /7 _1/2 (-// _1 /7^ 1 +X 2 - ^- 2 ))T n (/i 1 /2/,-l/2 ( _ ;U -l^ l +y2 _ ^-2)) e /6(x 1 - W ) ^ 

Plugging ^ = /i/)^ 1 (|(x 2 + y 2 ) + C — v ^ 2 ) one can rewrite (3.8) as 
(3.9) e(x, y, 0) = (2tt)- i / ,3/2 /) -3/2 e / M /,-i(i(x 2+y2 )-^- 2 )(xi-yi) x 



E 



x 



n J(2n+l)fih<W+v(x 2 +y 2 )+2vC-v2 t j,-2 

T n (/i 1/2 /)- 1/2 (^(x 2 - y 2 ) - C))T n (// 2 /r^(i( X2 _ y 2 ) + ^e-^-^-n) d( 

where factor e lfMh 1 ^H X2+ y 2 )~ v ^ 2 )( x i~w) cancels with the adjoint factor coming from e(y,x, 0) 
in the final calculation of / = / 2 . Note that W + v(x 2 + y 2 ) is potential V calculated at 

\{X2+Y2). 

In particular as v = (degenerate case) one gets 
(3.10) e(x,y,0) = e^ w (x,y, 0) * f (2,r)-y/ 2 r 3 <'M i »"'~ I < Xi+ >' ! ><"->'>> 



E 



x 



n J(2n+l)fih<W 



T n ^ l/2 h-^{\{ X2 - y 2 ) - C))T„(/, 1 / 2 /7- 1 / 2 (i(x 2 - y 2 ) + C )) e^*" 1 ^) c*C- 

Calculating the trace of this kernel we get (27r) _1 • #{n : (2n + < 1/1/} which is well 
known Magnetic Weyl expression. 

Note that (3.10) (modulo exponential-quadratic factors) depends only on |x — y| and 
one can rewrite it as 



x 

— J(2n+l)^h<W 



(3.11) en*,y,0) = (27r)-y /2 /r 3/2 V / 

n J (2n^ 

y n {^ 2 h-^{\\x - y| - cYn^ 1/2 h^ 2 {\\x - y| + C)) d( 
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or equivalently 

(3.12) CV.y, 0) = (27r)-V 3/2 ^ 3/2 e^'^ 1(x2+y2)(x44 - yi) x 

„ J(2n+l)iih<W ^ ' ^ ' 

3.2 Calculations with Oscillatory Integrals 

It immediately follows from [Ivrl] that as \ t\ < e/i propagator u is given by 

(3.13) u(x,y, t) = TlUTl = pit^h)- 2 J f e'^^^-^^-s^^W^) x 

E «/(*. 04(y. /. CK'-''/^'' E T n (x[^ 2 h-^)T4yy/ 2 h-^)u nn/ (x^ y' 2 , t) d^drj 

1,1' n,n' 

where T is /i _1 /?-FIO, 

(3.14) Fv{x) = (27r/i- 1 /))- 1 f e '^- 1 (S(x, 5 ,M- 1 )-(x',?» ^ a/ ( x _ x '_ £)„( x ') dx > d ^ 

which reduces A to its canonical form A = T*AT with 

(3.15) A = u{x 2 , f i- 1 hD 2 ){h 2 D 2 + /i 2 x 2 ) - W{x 2 , /i- 1 /)D 2 )+ 

E a mW (x 2 , ^hD 2 )(h 2 D 2 + ^r^-^-'h 1 

m+k+l>2 

with cj = FoM/, UV = V o \|/ with some map V : T * R -> M 2 . 
Really, decomposing 1/ (propagator of .A) into sum 

(3.16) U(x',y',t) = E ^(x2-y2-t)T n (x^ 1 /2/ 7 -V2 )Tn/(y / // i/2^-i/2 ) 
one gets (3.13). 

On the other hand, since A is a diagonal matrix in the basis of T n (x{/i 1 / 2 /? -1 / 2 ) I conclude 
that the non-diagonal elements u nn < are negligible on the time interval I are interested while 
diagonal elements are Schwartz kernels of propagators e lh Ant of 1-dimensional operators 

(3.17) A n = u{x 2 , f i- 1 hD 2 )({2n + l) f ih) - W(x 2 , /i' 1 hD 2 )+ 

E a mW (x 2 , /U - 1 /)D 2 )((2n + l)ph) m n 2 - 2m - 2k -'h' 

m+k+l>2 
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with the principal parts 

(3.18) A° n = u{x 2 , fJ ,- 1 hD 2 )({2n + l) f xh) - W{x 2 , ^ hD 2 )+ 

3 mkl {x 2 ,^ 1 hD 2 )((2n + l) f ih) m l 2 



2-2m-2/c 



m+k>2 

-1 A <J 



Rescaling t — > t' = fj, 1 t one gets e lflh Ant with \t'\ < e which are standard /i ^-FIOs; 
more precisely 

(3.19) u nn (x' 2 ,y' 2 , t) = {2n fi-'h)- 1 J e '^(M4.v^- 2 )+My^- 2 )t') x 

Pi{n, x' 2 ,y 2 ,t', /i- 2 )/^'/)-' | tl=fi _ lt d( 

i 

where 

(3.20) A° n (x 2 ,d 4 <f )n ) = A n (y,(), 

(3.21) <p n = 0, d x ^(f) n = ( asx' 2 = y 2 ,. 

So, one needs to plug (3.13) with sum over n = n' with u nn defined by (3.19)-(3.21) into 
Tauberian expression (1.4) for e(x,y,0) and then to plug the result into calculation of X. 
One can prove easily that skipping terms with / > 1 or /' > 1 in (3.19) and (3.13) results 
in 0(/i -1 /?~ 1 ~ K ) error in I. Also, with the same error one can replace T = e/x by T = oo 
after these substitutions are made leading to the following analogues of (3.13) and (3.19): 

(3.22) e'(x, y , 0) = (In^h)- 2 J j e^" (s^^Mx'.O-sty.^-^^o) x 

a(x, x', 0«(y, y\ j £ T^y^/i-^T^^/i-^e^, y 2 , t) 

with 

(3.23) e f ,(x^y^t) = (27r Ai - 1 /))- 1 / ^"^(x^c^^ x 2 , y 2 , 0, £ //~ 2 ) <*C 

7{^(x^,^,C)<0} 

where a = ao and (5 = (5q. 
So, I arrive to 

Proposition 3.1. Let conditions (1.1) — (1-3) be fulfilled. Then replacing in I e(x,y,0) by 
e'(x,y,0) defined by (3.22) — (3.23) we make an error 0(fi~ 1 h~ 1 ~ K ). 

However even if one can make more descriptive construction defining all the phases and 
amplitudes in geometric terms 9 -*, I prefer to make a less sharp but more explicit calculations. 



9 ) The reader who is more geometrically savvy than me can do this. 



3 CALCULATIONS: STRONG MAGNETIC FIELD 



30 



3.3 Calculations in Outer Zone 

3.3.1 Consider outer zone {\x — y\ > 71} where 71 > Co/i -1 will be defined later. It 
follows immediately from proposition 2.1 that 

(3.24) The contribution of zone {|x — y| > 71} with 7 > Co/i" 1 to the asymptotics does 
not exceed 

(3.25) Cfi- 1 h- 1 1 - 1 - K 

and therefore does not exceed C/i~ 1 /?~ 1 ~ K as 7 = 72 = h n ^ K+1 \ However I want a better 
estimate. 

Let us consider first uj(x) = 1 and ipj = ipj n which are 7-admissible functions with 
7-disjoint supports and 7 > Cofi" 1 . Let us apply formula (2.13) [Ivrl2] for X 

(3.26) /T 2 Tr / F tl ^ h -i Tl ^ h -i^ (xr(ti)XT(t 2 )^i, 7 ^2, 7 ,t 1 U{h + t 2 )j) dr 

^(ri,r 2 )eIR-' 2 V 7 

with ij) t = U(t)ipU(—t). Originally this formula is with However due to propagation 
for magnetic Schrodinger operator one can replace Xe^(0 by Xr(t) where T = /17 and x(0 
is supported in (— 2 Co, — e ) U (e , 2 Co) and equals 1 on (—Co, — 2eo) U (2eo, Co). 
After this substitution can rewrite (3.36) as 



T1=T2=0 



T = 



(3.27) -T 2 F tl ^ h -i Tlit2 ^ h -i T2 (^XT{ti)x T {t2)^i,^2, 1 ,t 1 U{t 1 + t 2 ) 

- (27T)- 1 7- 2 F t ^- lT (fx T {\t + s)xr{\t - s)Vi, 7 V> 2 , 7 i t+s *) U(t) 



r=0 



with x(0 = x(t)/t. 

First, let us calculate L/(t)^ 7 t/(— t). To do this let us go to the canonical form of 
operator /4 (see [Ivrl], Chapter 6): 



(3.28) ^= Z*Z + l^/(x 1 ,/i- 1 / ? D 1 )+ ^ b mn (x 1 ,/i- 1 /)D 1 )(Z*Z) m /i 

m+n+/>l 

with 

(3.29) Z = hD 2 + />x 2 , Z* = /?D 2 - />x 2 . 



— 2m— 2n— / r / 
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Then ip 1 (x) is transformed into if)L{x, \i 1 hD) which can be decomposed as 7 > \i 1+6 into 
an asymptotic sum 

(3.30) ^(x,yr x hD)~ Y. f,^i.^ 1 ''Di)(( ai ( ta T 

where ( = £2 + 1x2, C = £2 — 1x2 are symbols of Z and Z* respectively and w here means 
yU" 1 ^-quantization. Then (C ai C ) w is a symmetric product of a\ copies of /i _1 Z and ct2 
copies of \i~ x Z* . 

One can see easily that 

(3.31) |V^ 7 , Q | = 0( 7 - |Q| - |/31 ) Va,(3, 
and moreover 

(3.32) \^PU ^/' ai - Q2 g Qi g tQ2 ^ 7 )ovj/| = 0( 7 1 - |Qhl/31 ) Va,f3 

a! 

where 9 = X — iY and 5^ = X + iY, X and V are real vector fields and actually important 
is only that 55^ = — \FA g where A g is a positive Laplacian associated with the metrics 
(F-V)- 

Since [Z*, Z] = 2/i/? and [A, Z] = 2fihZ, [A, Z*} = -2fihZ* 

(3.33) U{t)ZU{-t) = e ^- 1 tz*z Ze -ih-^tz*z = e 2i>t z< 

U{t)Z*U{-t) = e *- 1 tz*z Z *e- ih ~ ltZ * z = e- 2i ^Z 

and I conclude that dropping all terms with a\ 7^ 02 and calculating an error one gets an 
extra yU~ 2 7 _1 factor in the expression T~ 1 gr{t) with 

(3.34) g T (t) = T- 1 J xr{\t + s) XT {\t - s)V>i, 7 V 2 , 7 ,| t+s <fe 

and in (3.27) itself; the latter would lead to 0(/z _1 /?~ 1 ~ K ) error in the final answer. 

Further, dropping terms with a.\ = 02 > n would lead to the extra factor 
the error. On the other hand, replacing Z*Z by F _1 1/1/ o v|/ brings even smaller error since 
I set r = 0. 

Under transformation of the operator to its canonical form (3.28), Qj is transformed 
into 

(3.35) Q'ix^'HD)- J2 &(*i./^ D i)(C ai C tQ T 

where one can prove easily using the same method as in the proof of (2.17), (2.18) of [Ivrl2], 
that 
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(3.36) Symbols g a satisfy the same inequalities (3. 31), (3. 32) as ip ltOL - 

Then to calculate (3.27) one can apply the standard approach of [Ivrl], Chapter 6 
resulting in the asymptotic decomposition 

/o Q7\ , ,-U-L-l-K \ ^ „ , ,-2n+/ -2n+m i / 

(3.37) ~ fi h 7 ^nmifJ' 7 " 

n,m,l 

where 

(3.38) kqoo = const / y?(x)6( W - 2 /i/j) c/x 

and other coefficients have the same form 10 ). 

This answer (3.38) has exactly magnitude /i _1 /?~ 1 7~ 1 ~ K as it should have; and since I 
am interested in the answer modulo 0(/x _1 /7 -1 7 _K ) I can skip any term with an extra factor 
7 (or lesser). In particular, in the above calculations I can replace »/) 27 i t+s by (ip a o (t> t ) w 

with 4> t = e M 1{Hw (because this error brings an extra factor /x _1 /77 _1 ); moreover, I can 
replace map <t» t by <t» xt = tHw ^ (because these maps coincide modulo 0(7)). 
It immediately implies 

Proposition 3.2. If in the outer zone {\x — y\ > n~ 1+5 } one replaces e(x,y,0) by (3.8) 

like expression for operator which in an appropriate coordinates has form (3.4) n ) the error 
would not exceed 0( / u _1 /? _1_K ). 



3.4 Successive Approximations 

Now I still want to cover zone {C /i _1 < |x — y\ < /i~ 1+<5 '}. For this I am going to compare 
operator in question and the model operator without using canonical form. 



3.4.1 I discuss successive approximation method which should be modified to the current 
problem. Consider first an abstract form. There is a perturbed operator A and unperturbed 
operator A and perturbation B = A — A Q with ||6|| = v. Let us consider U(t) = e lb ltA 
and U {t) = e ih ~ ltA °; then 

U{t) = U {t) + ih- 1 [ U Q {h)BU{t - ti) dt x 



10 ^ Actually they are sums of such terms with 6 replaced by 5^ which is not essential due to (1.3); I 
remind that W — V/F. 
X1 ' I leave to the reader to reach such expression in the arbitrary coordinates. 
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and iterating I get 

(3.39) U(t) = J2' kh ~ k I Uo{h)BU Q {t 2 )B ■ ■ ■ U Q (t k )BU Q (t - h t k ) dt, ... dt k + 

n h~ n f Uo(t l )BU Q (t 2 )B---U Q (t n )BU{t-t 1 t K )d tl ...dt K 

JA K (t) 



0<k<K 
i' 



'A K (t) 

where 

A k (t) = {t= (t! t k ) : ht- 1 > tfct" 1 > 0, (ti + • • • + t^t' 1 < 1}, 

term with k = is Uo(t) and the last term is negligible as \ t\ < T, 

(3.40) Tu < h 1+s 

and K is large enough. 

Let us consider term with 1 < k < n; one can rewrite it as 

(3.41) i k h~ k [ B tl B tl+t2 ■ ■ ■ B tl+ ... +tk dh ... dt k x U {t) 
with B 5 = U(s)BU{-s). Rewriting 

(3.42) B = J2 Ba,pZ a Z*?, [B^, Z] = [B a , p , Z*\ = 

a, (3 

I arrive to 

(3.43) B 5 = B a ^{fi' 1 s)e i ^ a '^ s Z a Z* 13 
and (3.41) becomes 

(3.44) i k h' k f V (B aiA ---B akA )(^ 1 t)e i ^dt 1 ...dt k x 



Z^Z*? 1 ■■■Z ak Z* /3k U (t) 

with at = («i a k ), (3 = (3 k ) and f(t) = J2i<j<k f i(E;</</c( tt /- A)- Taking Taylor 

decomposition of (B au/3l ■ ■ ■ B aki p k } (/U~ 1 t) and calculating this integral I get 



(3.45) i k h- k (s a pB ai ^---B ak i,(iM t ^+ E R^s^'-'t'e^ 



t'e' 1 "" | x 

(TeOT(«-/3),/<fc,s>0 

Z ai Z */3i ...Z ak Z*P k U Q (t) 



3 CALCULATIONS: STRONG MAGNETIC FIELD 



34 



where y\(ct — (3) is a finite subset of 7L k . 

Plugging into (1.4) I get a term in successive approximations of ej 

(3.46) ,*/r* 7-" 1 J Xr(t) (6^.^(0) ...S atlA (0)^+ 

5^ Rk a /3ais^ k+ '~ s t l e i ^ x Z^Z*^ 1 • • • Z ak Z*P k U Q (t) dt 

ae%a-f3)Cl k ,l<k,s>0 

which is equal to 

(3.47) / fc /7- fc 7- i yy*xr(0(6 a/3 e QlA (o)...e Q , 1A (o)^+ 

}^ Rkc&Tisl*~ k+ '~ s t'e i,MTt> ) x Z^Z*' 31 • • • z°"'Z*^e iTh ~ lt dtd T E {r). 

ae < n(a-/3)Cl k ,l<k,s>0 

Note that 

JJ XT^t'e^- 1 ^ dtd T E (r) = T l+1 J ^rTTT 1 ) d T E (r - a^h) 

with = t'x(t) = /t /_1 x(t). 

Therefore the trace norm of the first term in (3.47) does not exceed 

(3.48) Ch-^T"- 1 H \\B ajPj \\; 

l<j<k 

where I used an already mentioned inequality \\E (r) — E(r')|| 1 < C/? _2 (/?7" _1 + \t — t'\). 

On the other hand, as /x < e/? _1 the trace norm of the second term in (3.47) does not 
exceed 

(3.49) C^+'h-^T 1 - 1 J] \\B ajPj \\. 

l<J<k 

Then the contribution of zone {\t\ x 7"} and the mentioned terms to the final answer does 
not exceed (3.48) and (3.49) respectively with an extra factor which is equal to h~ K as 
T < C and to fi K T~ K as T > C . More precisely, contribution of zone { |x — y | > C h} is 
trivial and {|x — y| < C h} follows arguments of of the proof of proposition 1.2. 

Note that in my settings v = fi~ m (with m = 2 if F = 1 and I compare with operator 
considered in subsection 3.1) and that ||B a;/3j || < CjjT max ( m .M+/3/l). 
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Then the first terms in (3.47) result in the final contribution to the error not exceeding 

(3.50) Ch- 1 - k - K T k - 1 fi- mk ; 
further summation with respect to T < Cq results in 

(3.51) Or^-^-^l + Sfcillog/il). 

On the other hand, as T > C max(l,/i/?| \ogh\) one should replace expression (3.50) by 

(3.52) Ch' 1 -* T k - 1 fi- mk x (fi- 1 T)- K x Ch' 1 ' 1 * T k-i-n^-mk+ K . 

then summation with respect to T (from T = Cq to T = fi s ) as k = 1 results in 
Ch~ 1 ~ k fi~ k+K < C/j,~ m h~ 1 ~ K and all other terms are smaller. 

Further, if one considers \aj\ + \(3j\ > m + 1 at least in one of the factors then estimate 
acquires an extra factor so the final estimate would be 

(3.53) Ch- 1 - K f i- mk - 1 (l + 5 fc i| log h\) . 

On the other hand one gets from the other terms with / > 1 the same estimate (3.51) 
and with terms with / = one gets the value as T reaches its lowest value e/z -1 i.e. one get 

(3.53) as k > 2 and (3.51) as k = 1 (but without logarithm term). 

There is one special case of e(h\ log /?| ) 1 < /x < h^ 1 when one needs to sum (3.50) to 
the upper bound Cjih\ log h\ rather than to C but it does not affect the term with k = 1 
and all other terms are smaller. 

So, as F = 1 error estimate (3.51) with k = 1 is achieved: 

(3.54) Cfi- m h- 2 - K \ \ogh\. 
Probably one can get rid off logarithmic factor. 

3.4.2 Let us consider F which is not identically 1. Without any loss of the generality one 
can assume that F(y) = 1. Then let us consider perturbation B — (3hD t , with (3 = (F _1 — 1), 
bringing F to that case. Note that if T > Co then 

= O^Tx h\ log h\ T^ 1 ) = O^-^llog^l) 

since \\hD t \\ on the interval in question could be brought to Ch\ \ogh\ due to logarithmic 
uncertainty principle. Then ||6|| T < h 5 . 
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Let consider the whole interval (—7", 7") and plug successive approximations to Taube- 
rian formula (1.4) directly. Then hD t could be dropped on xr and I get 

(3.55) CT- 1 [ (F t ^- lrXr (t) f U Q (t - t')(3Uo(t') dt'dr 

J oo JO 

in the first successive perturbation term and the other look similarly. Then using approach 
of the previous subsection one can estimate k-th. term by 

(3.56) CT-H- 1 x { f i- 1 T) k - K < C f i~ 1 h- 1 - K 
and one should not worry about it. 

3.4.3 So far I assumed that fi < eh' 1 . As e < \xh < 1 1 can reduce to canonical form of the 
model equation and consider non-diagonal terms which are 0(/i _1 ) as perturbations which 
can be excluded, since ellipticity locally is broken for no more than 1 number "n" , leading 
to the diagonal perturbations 0(/i~ 2 //i/?) = 0([i~ 3 h~ 1 ) which in turn leads to 0(h~ 2 ~ K x 
fi^h^ 1 ) = 0(fi~ 3 h~ 3 ~ K ) error in the final answer; but in this case it is 0( / u _1 /? _1_K ). 

3.4.4 So, let us summarize. 

Proposition 3.3. (i) One can replace in (1.4) U(t) by Uo(t) with the final error not ex- 
ceeding C fi~ m h~ 2 ~ K ~ s . In particular, as fx > /) _1 /( m_1 ) -5 this error is 0(/x _1 /) _1_K ) and as 
jj, > /7-V™-<5 this error is 0(/)" 1 " K ). 

(ii) On the other hand, adding the second term of successive approximations one gets a 
final error CLi~ 2m h~ 3 ~ K ~ 5 . In particular, as ii > /)~ 2 /( 2m ~ 1 )- <5 this error is 0(/i _1 /? _1_K ). 

4 Calculations: Intermediate Magnetic Field 

Now I am going to combine methods of two previous sections. Namely, while in section 3 I 
did not consider points as distinguishable unless \ t\ > C , I will do it here using arguments 
of section 2. As a result in addition to threshold 7"* = min(z/~ 1 /) 1+<5 , C ) appearing in section 
3 another thresholds 7~ # = T+(/i, h, 0) and T = 7"_(/i, h, <fi) appear as J* 1 x ^ and x J T 
pairs are considered and for \ t\ > 7~# distinguishability arguments are used. This is helpful 
as 7"* > 7"*(0) at least for some q and 0. 
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4.1 Preliminary Analysis 

The previous section construction works as I am able to run successive approximations with 
T = C h~ 5 i.e. as v < h 1+s . As v = fi~ m this means /i > /j- 1 /" 1 -- 5 . Now I am going to 
investigate both cases 

(4.1) lr2 fi > h- l ' m -\ fi < h- 1 '" 1 -' 6 . 

As 7"* > h~ 5 is fulfilled I estimated the contribution of {\t\ > T*} by Cji^ l h~ 1 ~ K in 
subsection 3.3 already. So, as 7"* > h~ s I can reset in what follows 7"* to h~ 6 . 

Further, as 7"* = h~ 5 , /i < e(h\ log and C < T < T* one can estimate the 
contribution of \t\ x T by (3.52) which sums with respect to T to its value as T = Cq i.e. 
to Q^-^/r 1 -*. So, 

So, 

(4.2) Under condition (4.1)i one can redefine T* = C ; then contribution of zone {\t\ > T*} 
to the remainder estimate does not exceed 

(4.3) Cfi K - mk h- 1 - k + C^h- 1 -*. 

I will use more sophisticated approach (weak magnetic field approach) as \ t\ < C . Note 
that all our analysis makes sense as 7"* > efi^ 1 only i.e. 

(4.4) /i > /,-V(-+i)-5 

4.2 Analysis of J ± x J ± Pairs 

4.2.1 At this subsection I consider only x J ± pairs without polar caps; x J T pairs 
and polar caps I consider later. Let us pick up T = T*. Plugging 0x7 under condition 
(4.1)2, one can see that distinguishability condition 07 > (fih\ \ogh\) 2 ^ 3 is fulfilled as 

(4.5) 7> T* = C(fih\\ogh\y/ 3 . 

Note that (4.5) always follows from (4.4) provided m > 3. Further as m = 2 (4.5) follows 
from (4.4) as /i > /r 2 / 5 -' 5 . 

On the other hand, as m > 4 unperturbed operator is too complicated to handle and I 
am interested in the case m < 3 only. 

Proposition 4.1. Let condition (4.4) be fulfilled and let either m > 3 or m = 2, ji > 
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Then contribution to the error of zone \ t\ x T < C and J ± x P pairs to the error does 
not exceed 

(4.6) Ch- 1 fi K (T- 1 - K + T- 2K + T- 2 8 Kl \\ogh\) as T > T* 

(4.7) Ch- l - k ^- mk {T k - 1 - K + T k - 2K + b Kl T k ~ 2 \ log 7~|) as T* < T < T* 

(4.8) Ch- 1 - k fi- mk (fi K T k - 1 - K + T k - 2 h- K {fih\ log h\) 2/3 ) as 7* < T < 7* 

(4.9) Ch-'-^^-^T^ 1 asT = fi~ 1 <T< 7" 3 * 



where 
(4.10) 



r f C(/i/)| log/)|) 2 / 3 as /i > e(/7| log/)|)- 2 / 5 , 

3 ~ | e^ 1 as /i < e(/?| log /?|) 2/5 ; 



m the latter case zone (4.9) disappears. 

Proof, (i) Note that as T > 7~ 2 all P x ^ pairs are distinguishable and one can apply 
estimate (2.17): an error does not exceed 

(4.11) CT-'h- 1 J (/i-V7T K ^0; 

plugging T = T* one gets (4.6). 

(ii) For T < T* one can apply successive approximation method and (4.6) acquires factor 
fi~ mk h~ k T k thus resulting in (4.7). 

(iii) For 7~ 3 * < T < 7~ 2 * there is no distitinguishability on the ends {4> < 4>t = C(fih\ log h) 2 ^ T~ 1 } 
and integral in (4.10) is taken from (fi T to 1, while contribution of the ends is estimated by 
Ch~ 1 ~ k ~ K T k ~ 1 (f) T thus resulting in (4.8). 

(iv) Finally as T < 7~ 3 * there is no distinguishability even in the center and <p T x 1. □ 

Remark 4.2. Note that for m = 2,3 I assumed \i > /rVfm+iM an d then case (4.9) may 
occur. 

For m = 2, /j _1 /3-<s < /j, < /?~ 2 / 5_<5 the previous order 7~ 3 * <C 7~ 2 * <C T* is replaced by 
7~ 3 * <C 7"* <C 7~ 2 * and the same arguments imply 

Proposition 4.3. Let m = 2, h' 1 / 3 - 6 < fi < h~ 2 / b - 5 . 

Then contribution to the error of zone \ t\ x T < C and P x P pairs to the error does 
not exceed (4.6) as T > 7~ 2; 

(4.12) Ch^^T- 1 ^ + 7"- 2 /r K (/i/)|log/)|) 2 / 3 ) as T* <T< T*, 

(4.8) as T* < T < T* and (4.9) as t < T < T* ; as /i < e(/>| log />|)- 2 / 5 ^one (4.9) 
cfo appears. 
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4.2.2 Now I need to sum estimates (4.6), (4.7) or (4.11), (4.8) and may be (4.9) with 
respect to T in the corresponding intervals and find their sum. Note that (4.9) almost 
always sums to its value as T hits its largest value T = 7~ 3 * which is exactly the second 
term in (4.8) as T hits its lowest value T = T^; the singular exception is k = 1 when the 
summation of (4.9) results in extra term 

(4.13) C lf i~ m h' 2 - K (\og(^ 2 h\\ogh\)) 



+ 



In the analysis below I assume that k ^ 1; otherwise some terms may acquire extra loga- 
rithmic factors; details are left to the reader. 

(i) Let either k = 1 or k = 2, k > 1. Then almost every term in (4.6), (4.7) or (4.11), (4.8) 
is maximized by T hitting its lowest value; the sole exception is the second term in (4.7) in 
the case k = 1, k, < 1/2 but then it is dominated by the first one. Then summation results 
in (4.8) calculated as T = 7~ 3 * i.e. 



(4.14) 



Or 1 -^-™ VI lo S h\) 2{k - 1)/3 as fj, > e(h\ log h\y 2 / 5 ; 

Ch- l - k ^- k - mk {f K + Lih- K (/ih\ log h\) 2/3 ) as n < e{h\ log h\)- 2/5 . 



(ii) Let k > max(2, 1 + k,,2k), k < 1 (equalities would bring some logarithmic factors). 
Then every term in (4.7)-(4.9) is maximized by T hitting its highest value; then as 7~* < 1 
the main contribution is delivered by the first term in (4.6) as T = T* or by the second 
term in (4.8) as T = T* and I get C^- m ~ mK+K h- 2 - K ~ s and Ch- l - k - K ^ mk (iih\ log /j|) 2/3 ; in 
the opposite case as /i > /j- 1 /™-*' the contribution of zone {\t\ > T*} is given by (4.3) and 
again it is larger than contribution of other zones with the possible exception of (4.8); so 
the final answer is 



(4.15) 



c ^- m - mK + Kh -2-K-5 + Ch-^^^-^ifihl log h\) 2/3 asfi< eh- 1/m - 5 '; 

Cfi K - mk h- 1 - k + Cfi' 1 h~ 1 - K + Ch^- k - K fi~ mk {fih\ log h\) 2/3 as /i > eh' 1/m ~ 5 ' . 



(iii) Let k = 3, k > 3/2. In this case in (4.7) the second term is dominant but as T = 7" 2 * 
it is still less than the second term in (4.8); so the final answer is given by (4.15) again. 
So, I arrive to 

Proposition 4.4. In frames of proposition 4-1 the contribution of x J T pairs to an error 
does not exceed (4.14) in the case (i) or (4.15) in the cases (ii)-(iii). 

In frames of proposition 4.3 the main contribution is delivered either by (4.12) or by 
(4.8). Contribution of (4.12) always is its value as T = T* = fi 2 h 1+s , namely 

(4.16) C^ 2 ~ K h- 2 - K - & + Cfi- A/3 h- 7/3 - K - 5 . 
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Note that value of (4.8) as T = T* is less than this, and its value as T = T 3 * is 

(4.17) Qi^/r^V^I log hl) 2 ^- 1 -^ 3 + Cir 2k h- x - k - K {nh\ log h\) 2 ^-^ 3 . 
So, I arrive to 

Proposition 4.5. In frames of proposition 4-3 the contribution of x J T pairs to an error 
does not exceed (4.16) + (4.17). 

4.3 Analysis of J ± x J T Pairs and Polar Caps 

4.3.1 Consider now x J T pairs. Polar caps we consider later. As 7"* > C I already 
estimated contribution of {\t\ > T*} to the error. In this case contribution of T < T* 
to the error does not exceed Cfi~ mk h~ 1 ~ k ~ K and no improvements is possible. Really, as 
analysis of subsubsection 2.2.3 as T% > C with 

(4.18) T* = C/i 3 h\ log h\ 

weak magnetic field approach does not work at all. However 7"* = fi m h 1+s < T% as m < 3. 
This estimate is larger than I got for J T x pairs previously and so 

(4.19) As n > h-^ m ~ s the (prefinal) answer is C/i^/r 1 -^. 

On the other hand as 7"* < C one can use results of subsubsection 2.2.3 for T* < T < C 
and one gets Ch^ 1 ^ K \ log h\ as // > e(h\ log Z?!)" 1 / 3 . And since weak magnetic field approach 
here really works only as T > T%, the contribution of interval [7"*, T%\ to the error is 
estimated by Ch~ 1 ~ K \ log T%/T*\ x C/i _1_K | \ogh\ and this estimate cannot be improved by 
our methods. 

Meanwhile contribution of J ± x J T pairs and T < T* to the error is smaller than 
Ch- l - K \ \ogh\ as k > 1. 

So I arrive to the same error estimate as before: 

(4.20) As /i- 1 /( m + 1 )-* < fi < h-^ m ~ s the contribution of J ± x J T pairs to the error does 
not exceed C/7 -1 ~ K | \ogh\. 

I remind that the case jx < (h\ log /?|) -1 / 3 should be addressed only as m > 3. 
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4.3.2 Let us consider polar caps {|sin0| < p = C(ph\ log Z?!) 1 / 2 }. As \t\ x T their 
contribution to an error does not exceed 



(4.21) Ch~ 1 ' K T- 1 x (p^h' 1 T) k x p x 

C/?- 1 ~>- m+1/2 /?- 1/2 | log h\- 1/2 (fi~ m h- 1 T) k -\ 



Summation of (4.21) with respect to T acquires an extra logarithmic factor as k = 1 
but /i-^ 1 / 2 /?- 1 / 2 << So 

(4.22) As /)~ 1 /( m + 1 )-< 5 < ^ the contribution of polar caps to the error does not exceed 
Ch~ 1 ~ K \ \ogh\ as well. 

4.3.3 So, as m = 2, 3 and / 7 - 1 /( m + 1 )+' 5 < ^ < /,-!/">-« the contr ibution of > x pairs 
and polar caps to an error does not exceed Ch~ 1 ~ K . This result is valid for k = 1 and is not 
improved as k > 2. I need to add this the contribution of J^ 1 x pairs arriving to 

Proposition 4.6. Let m = 2, 3 and < ^ < 77^ 

(a) In frames of proposition an error does not exceed (4.14) + Ch~ 1 ~ K in the case (i) or 
(4.15) + C/?~ 1_K in the cases (ii)-(iii) of sub sub section 4-2.2; 

(b) In frames of proposition 4-3 an error does not exceed (4.16) + (4.17) + Ch^ 1 ^ K . 
I remind that as /i > /j- 1 /" 7 ^ 5 an error does not exceed Cfi~ mk h~ 1 ~ k ~ K . 

4.4 More Calculations 

At this moment all calculations are done in the case m = 2, k = (so actually estimate is 
achieved with k = 1; see subsection 3.1). I am going to consider m = 2, k > 1 leaving to 
the reader to consider case m = 3, k = 0: in this case unperturbed operator is quadratic. 

So, in order to get above estimates one actually needs to calculate the unperturbed 
expression and the next term in the successive approximations. 

To exploit perturbation, let us note that perturbation to be considered is 



Here one needs to consider only T < p. In this case 



p- m h- 1 T< p- m h- 1 {ph\ \ogh\) 



/ 2 ^{p l - 2m h- l \\ogh\f' 2 <h 5 



as p > /,-V(2 m +i)-5 



(4.23) 
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and m = 3 in our case because an extra factor (x — y) results in an extra factor /i 1 . 
The k-th perturbation term in the decomposition of E(0) is equal to 

(4.24) -h' 1 -" f F t ^ h -Jx T (t)x 

J —CO 

U (t - tJBUofa - t 2 )BU (t 2 - t 3 ) • • • U (t k ) dh ■ ■ ■ dt k ^j dr = 

- h' 1 -" J° F t ^ h -i T (xT(t)U (t) B tl B t2 ■ ■ ■ B tk t/ti • • • dt k ) dr 

where A t = {(ti, ... , t k ), < ±t k < rkt k -\ < ••• < ±ti < ±t} depending on the sign of t 
and = Uo{-t')BUo{t'). 

To analyze (4.24) let us consider U (t)(xj — y,)l/ ( — t). Note that 

(4.25) xj = x'j + Xjii' l Z + xj/^Z* 
with slow evolution of the components 

(4.26) U {t)x'jU {-t) = x'j + fi^vjt + 0{fi- 2 ), 

and U {t)ZU {-t) = e 2i ^Z, U {t)Z* U {-t) = e'^Z (see (3.33)). 
Let us decompose 

(x- y r= v- m p^(x'-yr(z*-z y r(z x -z y y^ 

|/3| + M=3 

where a, [3,a £ Z +2 . 

As (x' — y') applies to 5(x — y) one gets 0; therefore modulo term which brings an extra 
factor /i -1 and leads to a final error 0([i~ l h~ l ~ K ) even as k = 1 

u (t)(x-yru (-t) = v- 3 £ p^(^x-z y ) /3l ( e - 2 '> t z;-z;) /32 (^r. 

I0I+M=3 

Opening parenthesis and integrating with respect to t one can skip with the the same 
final error 0(/i _1 /? _1_K ) all the terms containing Z x or Z* with the exception of the case 
Pi = P2 = 1; so only the only important term is 

l*(t)(x - y) a Uo(-t) = f^- 3 £ (-l) l/,| ^ZA^A( vt )-+ 

|/3| + k|=3,/3^(l,l) 

kl=i 
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One can see easily that the order of Z y , Z* in the product does not matter. Then 

B tl ■ ■ ■ B tk = ^ k p^,nZ^z^{z;z;yv^ 

|/3|+H+2n=3/c, \u\=\a\ 

with t = (ti, ... , t k ) and v e Z +k . Then integration over A t results in 

B tl . . . B tk = ^ Yl Pk ,,,nZ^{z:z:ye k -\^ 

\/3\+2n<3k 

and then (4.24) is equivalent to 

//- 3fc ^H^-> M , CT ZfZf(Z;z;)"9f-l^- 2 "Eo(x,y,r)| r=0 

\/3\+2n<3k 

and the total perturbation term in E(x,y,0) is 

(4.27) £ £ ^^^'^^ZfZflZ^^^-I^^Eolx.y.r)!^ 

l<y<fc-l |/3|+2n<3/ 

and total perturbation term in the final answer is 

(4.28) Yl fi~ 3J h i ~ W ~ 2 "^ 

l<J<k-l \/3\+2n<3j 

J j Pj ^x,y)Z^Z^{Z;Z;rd 2 J^\- 2n E {x,y,T)\ r=Q x E (y, x.O) dxdy. 
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